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Plasma  Current  Sheet  Development  in  a 
Clamped  Discharge  Thruster 

C.L.  Dailey  and  R.H.  Lovberg 


1  Introduction 

This  report  was  prepared  by  the  Proptilsion  and  Power  Laboratory  of  TRW’s  Applied 
Technology  Division.  It  describes  work  performed  for  the  Air  Force  Office  of  Scientific 
Research  under  Contract  No.  F49620-87-C-0059,  to  investigate  the  ionization  and  current 
sheet  development  in  an  initially  im-preionized  gas  and  to  evaluate  the  effect  of  clzonped 
discharge  operation  on  thruster  efficiency. 

The  Pulsed  Inductive  Thruster  (PIT)  used  in  this  investigation  produces  the  plasma 
ciirrent  inductively  rather  than  by  an  electrode  discharge.  It  has  been  investigated  pre¬ 
viously  [1]  as  a  candidate  for  Electric  Orbit  IVansfer  Vehicle  (EOTV)  propulsion.  In  the 
course  of  technology  development,  supported  by  the  Air  Force  Astronautics  Laboratory, 
the  question  arose  whether  a  clamped  discharge  would  be  more  efficient  than  the  ringing 
discharge  mode  that  had  been  used  for  all  of  the  prior  work.  Surprisingly,  the  result  was 
an  unconvincingly  firm  negative.  One  of  the  changes  made  in  converting  the  thruster  for 
clamped  discharge  operation  was  to  mount  it  coaxially  with  the  vacuiun  chamber  with  a 
relatively  small  separation  between  the  coil  and  the  chamber  wall.  Theoretical  calculations 
of  the  wall  interference  effect  indicated  that  this  could  have  accounted  for  the  poor  perfor¬ 
mance.  To  eliminate  this  possibility,  the  one  meter  diameter  thruster  coil  was  replaced  by 
a  67  cm.  coil  for  the  present  study.  This  also  resulted  in  a  50%  increase  in  the  azimuthal 
electric  field,  which  increased  the  range  available  for  the  present  research. 

The  results  of  the  experimental  work,  and  theoretical  modeling  performed  in  order  to 
help  understand  these  resxilts,  are  presented  in  this  report. 

2  Summary 

A  new  thruster,  designated  MkIV,  has  been  constructed  and  tested.  It  is  characterized  by 

1)  a  snudler  diameter  coil  than  its  predecessors,  i.e.,  67  cm,  as  compared  to  1  meter,  and 

2)  a  provision  for  diode  clamps  in  its  capadtor  bank,  which  prevent  reversal  of  the  driving 
current,  and  thereby  restilt  in  much  increased  capacitor  life.  Measurements  of  and  effi¬ 
ciency  have  been  made  on  an  impulse  balance,  for  both  clamped  and  undamped  operation. 
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and  a  detailed  survey  of  the  electromagnetic  field  structure  and  plasma  properties  in  the 
accelerating  region  has  been  made  for  the  undamped  mode. 

Thrust  balance  measurements  showed  the  effidency  of  this  thruster  to  be  significantly 
lower  than  Mkl,  and  no  better  than  Mklll.  Undamped  (rinpng-mode)  operation  was 
significantly  more  effident  than  the  damped  mode. 

Probe  measurements  of  the  electric  and  magnetic  field  and  current  density  distributions 
allowed  a  determination  of  the  ionization  levd  of  the  plasma,  and  inference  of  the  effidency 
of  the  induced  plasma  current  layer  as  an  accelerating  piston.  It  was  found  that  the  MklV 
propellant  was  very  weakly  ionized,  and  only  partially  accelerated  by  the  first  current  sheet. 
This  is  in  contradistinction  to  Mkl,  which  displayed  tight  current  sheets,  and  complete 
ionization. 

A  ntunerical  simtilation  of  the  plasma  breakdown  and  current  sheet  formation  process 
allowed  a  comparison  between  the  Mkl  and  MkIV  thrusters,  based  on  such  parametric 
differences  as  operating  voltage,  propellant  density,  and  parasitic  circuit  inductance.  It 
was  found  possible  to  achieve  quite  good  agreement  between  predictions  of  the  code  and 
observed  behavior  of  the  two  thrusters.  It  could  be  concluded  that  the  inefficiency  of  MkIV 
was  due  almost  entirely  to  its  parasitic  inductance  being  too  high.  The  effect  is  that  if 
the  current  rise  rate  is  not  sufficiently  high,  electron  temperature  is  kept  depressed  by 
radiation,  and  the  ionization  rate,  exponentially  dependent  on  temperature,  is  too  low  to 
provide  a  highly  ionized  plasma  in  the  time  of  the  current  rise. 

This  deficiency  in  the  design  of  MkIV  is  an  easily  correctible  one;  redesign  of  the 
connecting  and  switching  circuitry,  as  well  as  the  use  of  multiple  parallel  low-inductance 
capacitors  should  bring  the  relative  parasitic  inductance  of  the  system  to  a  value  below 
that  of  Mkl. 

3  Background 

The  Pulsed  Inductive  Thruster  is  one  of  three  electric  thrusters,  capable  of  performing 
federal  and  military  spacecraft  missions,  that  use  electric  or  magnetic  forces  to  accelerate 
the  propellant.  The  others  are  the  Ion  Engine  and  the  MPD  engine,  the  coaxial  counterpart 
of  the  inductive  thruster  which  operates  in  the  steady  mode. 

The  PIT  is  inherently  a  pulse  mode  device.  It  differs  firom  other  (non-thermal)  thrusters 
in  that  it  has  no  electrodes,  and  from  other  inductive  accelerators  in  the  use  of  a  flat, 
spiral  coil  that  forces  constant  radial  current  density  in  the  coil.  While  prior  technology 
work  showed  the  PIT  to  be  capable  of  achieving  efficiencies  high  enough  for  practical 
applications,  the  poor  effidency  attained  during  clamped-mode  experiments  showed  also 
that  the  conditions  necessary  to  achieve  high  effidency  were  not  well  tmderstood. 

The  1  m  diameter  ringing  mode  thrusto’,  Mkl,  used  for  the  first  thrust  balance  investi¬ 
gations,  gave  the  highest  effidency  and  by  far  the  most  highly  diamagnetic  plasma  current 
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sheets  of  any  device  tested  [1].  This  thruster  was  not  operated  in  the  clamped  mode.  A 
revised  version  of  this  first  thnister  (Mklll)  was  tested  in  the  clamped  mode,  but  not  in  the 
ringing  mode.  Its  efficiency  was  found  to  be  far  inferior  to  that  of  Mkl.  It  was  concluded, 
however,  that  the  clamping  itself  was  not  responsible  for  the  change,  since  during  the  first 
quarter-cycle  of  current,  when  most  of  the  momentum  transfer  takes  place,  the  clamping 
elements  had  not  gone  into  conduction. 

The  present  program  is  a  continuation  of  this  work,  in  which  a  new  coil,  designed  to 
improve  the  strength  of  the  ionizing  electric  field,  was  incorporated  into  the  thruster.  A 
more  extensive  diagnostic  and  analytical  effort  has  been  undertaken,  in  an  effort  to  better 
understand  the  imderlying  physics  of  the  plasma  formation  and  acceleration  process,  and 
thereby  to  have  a  better  basis  for  efficient  thruster  design. 

3.1  Advantages  of  the  Concept 

FVeedom  from  electrode  erosion,  a  unique  feature  of  this  thruster,  ensures  very  low  space¬ 
craft  contamination  and  makes  it  possible  to  design  a  thruster  of  unlimited  life.  Because 
of  the  large  coil  size  relative  to  the  MPD  thruster,  its  radiant  power  loading  is  significantly 
less;  and  because  of  its  small  thruster  area  relative  to  the  Ion  Engine  it  is  capable  of  han¬ 
dling  more  power  for  a  given  engine  size.  Consequently,  it  is  a  natural  choice  for  a  high 
power  EOTV.  The  pulse  mode  of  this  thruster  makes  it  well  suited  for  missions  requiring 
fine  impulse  bit  control,  and  for  applications  to  spinning  satellites  in  which  the  impulse 
bits  must  be  sychronized  with  the  rotation. 

Another  advantage  of  pulsed  mode  operation  is  that  the  performance  of  the  thruster 
is  characterized  completely  by  the  properties  of  a  single  pulse.  Thrust  at  a  given  1,^  can 
be  varied  from  an  arbitrarily  small  value  to  the  maximum  available  for  a  given  power 
by  simply  increasing  the  pulse  rate.  Pulsed  mode  capability  is  also  a  great  advantage  in 
experimental  work.  The  sensitivity  to  background  pressure,  that  plagues  D.C.  thrusters, 
is  entirely  absent  for  the  inductive  thruster  which  produces  its  characteristic  impxilse  bit 
before  gas  rettuning  from  the  wall  has  had  time  to  reenter  the  discharge.  In  this  way,  the 
thruster  operates  in  an  effectively  perfect  vacuum,  even  with  slow  facility  pumping  speed. 

3.2  Problems  of  the  Concept 
3.2.1  Propellant  Iiyection 

Proper  placement  of  the  injected  gas  is  difficult  to  achieve.  For  high  efficiency,  the  gas 
layer  must  be  located  dose  to  the  coil,  and  a  high  fraction  of  the  gas  injected  must  be  in 
that  closely  coupled  layer.  Work  done  to  date  in  this  area  has  been  adequate  to  permit 
studies  of  ionization  and  current  sheet  btiildup,  and  thrust  balance  measurement  of  I«p  and 
effidency,  but  for  a  practical  thruster  more  effort  will  be  required  to  improve  the  radial 
uniformity  and  tightness  of  coupling  between  the  mass  layer  and  the  coil. 
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3.2.2  Propellant  Valve 

The  gas  injection  valve  is  a  magnetically  driven,  mechanical  device  that  must  operate  in 
a  fraction  of  a  millisecond  for  an  essentially  unlimited  cycle  life.  In  the  earlier  program 
[1]  a  life  of  a  few  million  shots  was  demonstrated  with  a  valve  design  that  used  a  thin, 
stainless  steel  diaphragm  brazed  to  a  cylindrical,  alumina  coil  form.  Repeated  frdlures 
were  encountered  due  to  fatigue  cracks  associated  with  stress  concentrations  induced  by 
the  brazing  procedure.  A  different  method  of  attachment  is  needed  that  does  not  rely  on 
brazing. 


3.2.3  Switching 

Experience  with  the  1-meter  diameter,  ringing  discharge  thruster  has  shown  that  a  highly 
diamagnetic  ciirrent  sheet  can  be  ionized  and  built  up  in  a  time  small  compared  to  the 
ourent  rise  time  of  the  first  half-cycle.  However,  it  is  essential  that  the  current  rise  time 
must  not  be  limited  by  the  characteristic  time  of  the  switch.  In  this  regard,  the  spark 
gap  is  ideally  suited  for  laboratory  work.  Its  switching  time  is  very  short  and  single  shot 
operation  is  entirely  satisfactory.  Conceivably,  spark  gap  switching  might  be  feasible  for  a 
practical  thruster;  but  a  potentially  more  attractive  concept  is  an  assembly  of  solid  state 
switches  known  as  Reverse  Biased  Diode  Thyristors  (RBDT).  No  work  toward  developing 
a  solid  state  thruster  switch  has  been  done  as  yet. 

3.2.4  Energy  Discharge  Capacitors 

This  has  been  regarded  as  a  difiicult  technology  area  in  the  past.  However,  current  state  of 
the  art  capacitors  are  now  capable  of  virtually  imlimited  discharge  life  at  an  energy  storage 
density  of  15  joules  per  pound.  Test  screening  at  several  million  discharges  is  effective  in 
improving  lifetime  reproducibility. 

3.2.5  Thruster  Mass 

For  high  power  mission  applications,  the  power  related  mass  of  an  EOTV  is  dominated  by 
that  of  the  power  source.  The  combined  mass  of  thruster  and  power  processor  is  about 
the  same  for  the  three  high  I«p  electric  thruster  candidates. 

4  Experimental  and  Theoretical  Background 
4.1  Experiments 

For  several  years,  TRW  has  conducted  experiments  with  pulsed  inductive  accelerator  coils 
ranging  in  diameter  from  0.1  m  to  1.0  m.  The  earliest  work,  supported  by  AFOSR, 
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iised  static  fill  gas  loading  with  coil  diameters  of  0.1,  0.15,  0.2,  and  0.25  m.  Diagnostics 
consisted  of  miniattire  plasma  probes  (  magnetic  and  electric  fields  and  current  density). 
Subsequent  work,  supported  by  AFRPL,  used  larger  coils  0.3  and  1.0  m  in  diameter  and 
thrust  balance  data  in  addition  to  the  plasma  probe  measurements.  Both  static  fill  and 
pulsed  injection  gas  loading  were  used  with  these  coils,  the  latter  being  necessary  for  the 
thrust  balance  measurements.  In  all  of  this  work,  the  accelerator  circuit  consisted  of  the 
coil,  an  energy  discharge  capacitor  bank,  and  a  spark  gap  switdi.  In  every  case,  regardless 
of  thruster  size,  a  highly  diamagnetic  current  sheet  formed  in  the  \m-preionized  gas  in 
about  200  nanoseconds,  a  time  small  compared  to  the  2  microsecond  current  rise  time  of 
the  accelerator. 

The  reason  for  going  to  coils  of  larger  diameter  was  that  the  efficiency  was  observed  to 
increase  dramatically  as  the  diameter  increased.  Efficiency  was  practically  undetectable 
for  the  0.1  and  0.15  m  diameter  coils.  For  the  larger  coils  it  was  still  very  low  but  rapidly 
improving  with  increasing  diameter.  The  results  were: 

Diameter  Igp  t} 

0.2  m  1200  s  0.05 
0.3  m  1470  s  0.18 

Table  1:  Itj,  and  efficiency  for  two  thrusters. 

A  theoretical  model,  anchored  to  these  data,  was  developed  to  indicate  the  efficiency 
improvement  that  might  be  expected  by  increasing  the  diameter  to  1.0  m,  the  largest  size 
that  seemed  reasonable  for  the  available  1.22  m  vacuum  chzunber.  The  promising  results 
of  this  analysis  [1]  led  to  the  construction  of  two  1.0  m  diameter  coils,  one  of  copper  tubing 
supported  by  TRW,  the  other  of  a  multistrand  flat  ribbon  similar  to  Litz  wire  supported 
by  AFRPL.  Thrust  data  were  measured  only  for  the  latter  coil. 

5  Thruster  Parameters 
5.1  Coil 

As  mentioned  earlier,  the  present  thruster  coil  is  similar  to  all  of  the  previous  piilsed 
inductive  accelerator  coils.  It  differs  only  in  diameter  which  is  67  cm.  The  smaller  diameter, 
compared  to  the  1  m  coil,  provides  a  greater  azimuthal  electric  field  for  a  given  capacitor 
voltage  and  eliminates  the  possibility  of  inductive  coupling  with  the  1.22  m  diameter 
vacuum  chamber  wall.  The  coil  coxisists  of  24  (me  turn  spirals  separated  15  degrees  from 
each  other  and  connected  in  parallel.  Current  is  returned  through  radial  conductors  on 
the  back  side,  separated  1  mm  from  the  spirals.  The  resulting  current  density  in  the  coil 
is  purely  azimuthal  and  constant  everywhere  except  over  the  outer  4  cm  radius  where  the 


strands  are  more  closely  spaced  to  increase  the  local  radial  magnetic  field.  This  technique 
was  used  with  the  previous  1  m  coil  to  compensate  for  the  loss  of  magnetic  pressure  that 
would  normally  occur  at  the  c(^  edge  due  to  fidd  frin^ng. 

The  coil  strands  are  placed  in  grooves  machined  in  plastic  plates  and  the  sandwich 
impregnated  with  epoxy.  This  leaves  the  coils  separated  1  mm  from  the  plasma.  Additional 
insulation  is  provided  by  a  3  mm  thick  sheet  of  glass.  The  resulting  4  mm  separation 
between  the  coil  and  the  plasma  represents  a  parasitic  inductance  that  is  small,  but  not 
negligible,  compared  to  the  vacuum  coil  inductance  of  about  520  nh.  While  this  separation 
is  much  smaller  than  the  8  cm  decoupling  distance  of  this  coil,  it  has  a  noticeable  effect 
because  of  the  high  magnetic  field  near  the  face  of  the  coil. 

The  spark  gap  configuration  of  this  thruster  is  also  different  from  that  of  the  1  m 
thruster.  It  uses  axial  rather  than  radial  ciurent  return  conductors  that  result  in  greater 
parasitic  inductance.  Other  differences  contributing  to  more  parasitic  inductance  are  the 
use  of  fewer  capacitors  and  a  smaller  number  of  coil  strands.  The  addition  of  diode  clamps 
inside  the  capacitors,  as  described  below,  also  added  some  parasitic  inductzmce. 

5.2  Capacitor  Bank 

The  capacitors  used  for  the  present  program  were  taken  from  the  earlier  1  m  thruster 
program.  The  original  1  m  thruster,  Mark  I,  had  a  single  spark  gap  and  four  5  fiF  capacitors 
for  a  total  bank  of  20  fjF.  The  second  version,  Mark  II,  used  12  smaller  capacitors  for  the 
same  capacitance,  each  with  its  own  spark  gap.  Toward  the  end  of  the  previous  program, 
12  of  these  capacitors  were  modified  for  clamped  discharge  operation  by  removing  some 
of  the  pads  to  make  room  for  the  string  of  diodes  used  for  the  clamp.  This  reduced  the 
capacitance  of  each  unit  to  1.34  /iF.  The  individual  spark  gaps  were  removed  and  the  12 
capacitors  discharged  in  parallel  through  the  single  spzu-k  gap  that  had  been  used  with  the 
Mark  I  thruster. 

This  modified  version  of  the  original  1  m  thruster,  Mark  III,  was  tested  at  the  end 
of  that  program.  A  few  thrust  balance  tests  were  run  up  to  20  kV  before  diode  failiires 
occured  in  two  of  the  units.  This  was  the  first  indication  of  very  poor  efficiency  associated 
with  clamped  operation.  It  was  also  the  only  thrust  balance  data  obtained  with  the  Mark 
III  thruster.  Consequently,  the  Mark  I  data  represents  the  only  thrust  balance  data  in 
the  ringing  mode  for  ai^  1  m  thruster.  This  leaves  open  the  possibility  that  the  poor 
performance  of  the  clamped  Mark  III  thruster  may  not  have  been  due  to  the  clamped 
mode,  but  possibly  to  some  change  in  the  thruster  configuration  that  was  the  same  for 
both  the  Mark  III  1  m  thruster  and  the  present  67  cm  thruster,  MkIV.  This  possibility  is 
discussed  further  below. 

Tests  were  run,  in  the  current  program,  \ising  the  10  remaining  1.34  ftF  clamped  ca¬ 
pacitors  at  voltages  up  to  18  kV,  while  6  additional  units  were  being  reworked  with  higher 
voltage  diodes  ciq>able  of  operating  up  to  25  kV.  These  high  voltage  clamps  were  physi- 
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cally  larger  than  the  original  ones,  requiring  further  reduction  of  the  capacitance  to  0.97 
fiF  each.  Because  of  the  high  cost  of  the  diodes,  only  6  capacitors  were  modified  for  25 
kV  clamped  operation.  After  testing  this  configuratimi  up  to  24  kV,  again  with  very  poor 
efficiency,  two  further  ciq>acitor  changes  were  made.  One  consisted  of  disconnecting  the 
diodes  in  the  ten  clamped  capacitors,  and  adding  two  of  the  original  nominal  1.67 fjF  units 
and  operating  the  MklV  thruster  in  the  ringing  mode  with  a  measured  capacitance  of 
17.4/iF. 

The  other  change  was  to  replace  the  pads  of  the  6  high  voltage  clamped  units  with  much 
higher  capacitance  pads.  By  going  to  a  lower  voltage,  all  film  pad  design,  the  capacitance 
of  each  unit  was  increased  to  4.5  fiF.  Because  of  a  failure  of  one  of  the  series  sections 
in  one  of  these  capacitors  during  acceptance  testing,  its  capacitance  increased  to  6.2  fiF, 
giving  a  total  bank  of  28.7  fuF.  This  large  change  in  capacitance  made  little  difference  to 
the  measiired  thruster  efficiency. 

5.3  Propellant  Increment 

The  propellant  mass  (argon)  injected  per  shot  ranged  from  1  to  10  milligrams.  The  pro¬ 
pellant  injection  valve  and  method  of  measuring  total  mass  and  mass  distribution,  both 
radially  and  axially  in  front  of  the  coil,  are  described  in  detail  in  Ref.[l]  .  Several  chiuiges 
were  made  to  the  gas  injection  nozzle  to  improve  the  distribution  over  the  smaller  coil. 
The  final  results  are  shown  in  Figure  1. 

6  Thrust  Balance  Performance  Data 

Figure  2  compares  thnist  balance  measurements  of  thruster  efficiency  and  for  the  1 
m  and  67  cm  thrusters  in  the  ringing  discharge  mode,  and  the  smaller  thruster  in  the 
clamped  mode  for  three  different  values  of  capacitance.  The  two  conclusions  to  be  drawn 
from  these  data  are  that  the  clamped  mode  is  much  less  efficient  than  the  ringing  mode 
for  the  smaller  thnister,  and  that  the  ringing  mode  for  the  smaller  thruster  is  much  less 
efficient  than  the  same  mode  for  the  larger  thruster. 

These  data  represent  the  total  thrust  balance  data  base  for  the  1  m  and  0.67  m  thrusters. 
Rirther  insight  into  their  performance,  obtained  from  an  analysis  of  probe  data  and  the 
theoretical  model,  is  described  later. 

7  Probe  Mapping 

Distributions  of  magnetic  field  (Br),  current  density  (7«),  and  axial  electric  field  (.E,)  were 
measured  by  means  of  appropriate  probes.  The  magnetic  probe  was  a  2  mm  diameter 
solenoid,  encased  in  a  6  mm  glass  jacket,  and  mounted  on  a  two-axis  carriage  that  allowed 
scanning  over  that  portion  of  the  r  —  z  plane  in  which  acceleration  takes  place.  Data  were 
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recorded  over  a  grid  of  60  points,  12  points  in  z  at  each  of  5  evenly  spaced  radii.  Z  stations 
were  spaced  at  5  mm  intervals  within  3  cm  of  the  coil  face,  and  at  1  cm  intervals  farther 
downstream.  Current  density  was  measured  with  a  miniature  Rogowski  coil  scanned  over 
the  same  spatial  grid. 

The  Br  probe  was  calibrated  by  direct  comparison  of  its  output  with  that  of  an  accu¬ 
rately  formed  1  cm’  standard  placed  in  the  same  pulsed  test  field.  Such  a  procedure  was 
not  possible  for  the  Rogowski  coil,  however,  since  its  senritivity  includes  its  effective  win¬ 
dow  area  as  a  factor;  this  area  will  generally  be  somewhat  different  from  the  geometrical 
area  of  its  opening.  Rather,  it  was  calibrated  against  an  axial  Br  gradient  measiued  in 
the  pulsed  plasma  at  a  position  near  the  coil  and  near  its  mid-radius,  where  the  dB^ldr 
contribution  to  is  negligible.  The  outputs  of  both  probes  were  integrated  by  a  passive 
R-C  integrator  of  120  /is  time  constant. 

The  axial  electric  space-charge  field  was  measured  with  a  differential  electric  probe. 
The  probe  electrodes  themselves  were  5  mm  lengths  #16  copper  wire,  aligned  along 
r,  and  separated  along  z  by  6.5  mm.  They  were  connected  to  the  primary  winding  of  a 
carefully  shielded  ferrite-core  pulse  transformer  whose  function  is  to  cancel  all  “common¬ 
mode”  components  and  leave  only  the  differential  signal. 

All  probe  data  were  recorded  by  a  digital  transient  recorder  of  8-bit  resolution  and  20 
MHz  sampling  rate.  Records  were  archived  on  magnetic  disk  media. 

8  Results;  Comparison  with  Mkl 

This  section  presents  distributions  of  the  density  and  electromagnetic  field  quantities  ob¬ 
tained  from  the  MklV  thruster.  Since,  as  has  been  mentioned,  the  eflBciency  of  this  thruster 
has  been  foimd  to  be  markedly  inferior  to  the  1-meter  Mkl  thruster,  data  from  that  device 
are  included  in  order  to  attempt  an  explanation  of  the  difference. 

For  these  diagnostic  nms,  the  bank  (17.4  /iF)  was  operated  without  diode  claunping. 
This  provides  a  twofold  advantage.  First,  the  dynamics  of  the  second  and  third  half-cycle 
plasma  currents  provide  an  indication  of  the  amount  of  gas  left  behind  by  the  first  current 
sheet  (where,  in  clamped  operation,  the  first  sheet  would  have  to  accelerate  the  entire 
gas  fill);  secondly,  an  informative  comparison  can  be  made  with  the  performance  of  the 
undamped  Mkl  thruster,  which  was  reas<niably  effident. 

The  data  shown  here  were  for  20  kV  on  the  bank,  and  an  ii^ected  propellant  mass  of 
3.6  X  10~'  kg  of  argon.  Since  the  physics  of  the  coupling  between  fields  and  propellant  was 
the  subject  of  this  run,  data  were  taken  only  at  a  radius  of  25  cm  (near  the  mid-radius 
of  the  c(^),  where  the  z-dei>endence  of  the  measured  parameters  could  be  assumed  to  be 
typical  of  the  coil  face  as  a  whole. 


8.1  B  Field  Distribution 


A  superposition  of  Br(f)  traces,  taken  at  axial  positions  z=l,..,6  cm,  is  shown  in  Figture 
3.  It  is  seen  first  that  all  the  traces  begin  to  rise  immediately  at  t=0;  this  is  indicative  of 
rapid  early  diffusion  of  the  field  through  the  plasma,  which  in  turn  indicates  poor  early 
conductivity.  The  negative  swing  of  the  traces,  which  be^ns  at  z=l  cm  at  about  8  ps,  is 
by  contrast  a  discretely  ordered  series  of  rises,  indicative  of  a  current  sheet  propagating 
outward  at  a  velocity  of  about  1.6  cm/pa.  With  the  magnetic  pressure  difference  indicated 
in  the  figure,  this  implies  an  average  plasma  density  encountered  by  the  second  current 
sheet  over  the  1-6  cm  range  of  about  5  x  10“m“®,  which  corresponds  to  approximately  20 
p  of  the  original  fill  pressure  of  approximately  600p. 

By  contrast,  a  corresponding  set  of  Br{i)  traces  from  Mkl  (Figure  4)  show  less  initial 
diffusion,  2ind  a  better  formed  first  current  sheet.  Subsequent  to  the  passage  of  the  current, 
the  traces  have  nearly  identical  time  behavior,  indicating  that  the  spatizil  field  distribution 
in  the  1-6  cm  zone  is  approximately  fiat;  this  can  only  occur  if  negligible  plasma  is  left 
behind  by  the  current  sheet. 

In  Figures  5  and  6,  Br(z)  is  shown  at  five  times  for  the  two  thrusters.  Here,  it  is 
also  evident  that  diffusion  of  the  accelerating  B  field  through  the  pl^lsma  is  much  more 
pronoimced  in  MkIV  than  in  Mkl. 

The  different  behavior  of  MkIV  and  Mkl  during  the  second  current  half-cycle  is  evident 
from  Figures  7  and  8.  In  the  8-11  ps  time  interval,  the  MkIV  Br{z)  distributions  exhibit  a 
strong  positive  slope,  indicative  of  substantial  mass  loading,  while  Briz)  for  Mkl  remains 
neeirly  fiat  as  it  swings  from  positive  to  negative  values. 

8.2  Azimuthal  E  Field 

A  premise  in  the  original  design  of  the  MkIV  coil  weis  that  its  reduced  diameter  (67  cm, 
as  compared  to  100  cm  for  Mkl)  would  incrcose  £$  by:  1)  decreasing  the  circumferential 
distance  over  which  the  terminal  voltage  would  be  applied,  and  2)  increeising  radial  spacing 
between  the  coil  and  the  vacuum  tank  wall,  so  that  the  effect  of  the  wall  as  a  shorted 
secondary  turn  near  the  coil  would  be  reduced. 

However,the  magnetic  probe  data  shown  in  the  previous  section  axe  suggestive  of  an 
actual  deficiency  in  Ee  in  MkIV,  since  early  conductivity  seems  lower  than  in  Mkl.  This 
hypothesis  can  be  tested  by  sui  integration  of  Br  over  z  and  differencing  such  integrals  at 
adjacent  timesteps  in  the  data. 

The  fiux  of  Br  per  imit  circumference  is 

*  =  f:B,(2.)Ax,  (1) 

isl 
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where  n  is  the  number  of  z  points  in  the  B  probe  scan.  Then, 


E8  = 


At 


(2) 


This  is  an  expression  of  Fwaday’s  law  that  holds  irrespective  of  the  ratio  of  diffusion  to 
convection  of  the  field  in  the  plasma. 

The  results  of  these  computations  are  shown  in  Figiires  9  and  10*.  It  can  be  seen 
that  for  MkIV,  the  initial  £$  is  approximately  1  kilovolt /meter,  an  unusually  low  value  for 
initiation  of  a  pulsed  discharge.  The  corresponding  initial  E$  for  Mkl  is  between  2  and  3 
kv/m,  the  larger  coil  diameter  notwithstanding.  In  each  of  these  curves,  the  subsequent 
rise  of  Ee  up  imtil  2  or  3  /is  is  due  to  the  addition  of  a  v  x  B  term  as  the  plasma  accelerates 
away  from  thr  coil. 


8.3  Plasma  Density 

If  separate  equations  of  motion  for  the  ions  and  electrons  in  a  plasma  are  written  and 
combined,  under  the  assumption  of  quasi-neutrality  £ind  scalar  pressure,  as  well  as  the 
negligibility  of  certain  terms  that  depend  upon  the  electron-ion  mass  ratio,  there  is  ob¬ 
tained,  upon  addition  of  the  two  equations,  the  MHD  equation  of  motion  for  the  plasma 
as  a  whole,  and  upon  subtraction,  a  relation  known  as  the  generalized  Ohm’s  law.  The 
derivation  is  carried  out  in  Spitzer’s  classic  monograph,  for  example[2]. 

A  form  of  Ohm’s  law  appropriate  for  these  experiments  is 

E-|-vxB-(l/n.e)[JxB-Vp.]-qJ  =  0  (3) 

where 

V  =  (rUeVe  +  miVi)/{mg  -H  m,) 

Pe  =  nfcTe 
T]  =  mJnt€^Tf 


where  r*  is  the  meem  collision  time  of  electrons  with  ions  or  neutrals. 

Since,  in  the  PIT,  v  is  purely  eixial  (along  z)  eind  J  is  along  0,  the  z-component  of 
Ohm’s  law  becomes 


E, - 

rige 


J$Bf.  — 


dpt 

dz 


=  0 


(4) 


In  the  event  that  the  J  x  B  force  term  is  dominant  over  Vpe,  this  reduces  to 


JeBr 

TlgC 


‘The  large  noise  at  early  times  in  these  results  is  due  to  the  differencing  of  adjacent  time  values  of  B, 
which,  for  the  first  microsecond  or  so  is  itself  a  noisy  signal. 
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which  then  enables  a  determination  of  the  electron  density  distribution  through  the  exper¬ 
imentally  measured  distributions  of  7#,  and  E^. 


(5) 


Previous  experiments  employing  these  combined  probe  analyses  have  usually  had  a 
domii  mt  Lorentz  force;  Equation  (5)  has  proved  adequate  to  produce  electron  density 
distributions  that  have  agreed  well  with  optical  refractivity  measurements  of  the  same 
plasma  [3]. 

Two  relational  criteria  must  be  met  by  B,,  and  in  order  for  Equation  (5)  to 
be  valid:  First,  since  E,  is  in  the  denominator,  its  passage  through  a  value  of  zero  must 
be  accompanied  by  a  simultaneous  passage  of  one  of  the  numerator  factors  through  zero; 
otherwise  an  infinite  density  would  be  imphed.  Secondly,  the  munerator  {JeBr)  amd  the 
denominator  (E*)  must  always  have  the  same  sign,  since  density  is  always  positive. 

Figure  11  is  a  superposition  of  raw  data  traces  of  7®,  and  E,,  taken  near  the  coil 
face,  at  z  =0.5  cm.  The  vertical  scale  corresponds,  in  each  case,  to  the  size  of  the  8-bit 
word  in  the  transient  recorder.  For  clarity,  E*  is  plotted  in  reverse  polairity.  It  can  be 
seen  here  that  the  zero  crossings  of  E,  are  nearly  simultaneous  with  crossings  of  one  or  the 
other  of  the  other  two  traces,  especially  for  times  earlier  them  10  ps.  The  stricture  on  sign 
agreement  is  also  well  obeyed.  Figure  12  is  a  calculation  of  ne(t)  from  these  data,  using 
Equation  (5).  The  short  transient  excursions  of  rig  below  zero  at  times  near  5,  13,  and  18 
ps  are  due  to  the  very  slight  residual  zero  crossing  disagreement.  The  three  distinct  peaks 
in  tig  correspond  to  the  formation  and  subsequent  sweeping  forward  of  the  three  successive 
current  sheets.  Their  relative  heights  are  indicative  of  the  fraction  of  the  original  propellant 
fill  left  behind  by  the  first  current  sheet. 

Figures  13  and  14  are  a  similair  set  taken  at  ^=1.0  cm.  The  peaks  in  n*  are  narrower 
and  higher  than  in  the  previous  set  because  the  pleisma  sheets  have  been  mewle  denser 
through  snowplowing  of  plasma  between  the  two  stations,  and  have  acquired  a  forwau-d 
convective  velocity  so  they  are  in  the  probe  position  for  a  shorter  time.  The  second  plasma 
sheet  is  particularly  thin  and  dense  because  it  is  compressed  between  the  two  opposing 
Br  fields  that  the  current  sheet  separates.  A  sharp  dip  may  be  seen  in  E,  at  the  time  of 
paissage  of  each  of  these  layers.  In  the  case  of  the  second  sheet,  the  dip  is  close  enough  to 
zero  to  introduce  a  large  probable  error  in  the  value  of  n*  inferred  from  it. 

In  Figires  15  amd  16,  sets  of  Bg,  7s,  amd  Eg  are  displayed  for  z=2.0  cm  and  3.0  cm. 
Here,  there  are  serious  departures  from  the  criteria  for  use  of  Equation  (5)  for  reduction 
of  the  data.  In  the  2  cm  case,  for  exaunple,  E,  reverses  sign  near  3  ps  and  11  ps,  where 
no  corresponding  crossing  of  either  B,  or  7s  occurs.  The  sharp  oscillatory  character  of 
the  fluctuation  in  E,  that  carries  it  through  zero  makes  it  clear  that  the  electron  pressire 
gradient  term  must  be  included.  If  a  sheet  of  significamt  electron  pressure  moves  forward 
in  z  past  the  probe,  the  signad  should  first  swing  negative  amd  then  positive  as  the  density 


11 


peak  passes.  The  clearest  instance  of  this  occurs  in  Figure  15  near  20  /is;  the  effect  is 
present  in  the  trace  near  4  and  10  pis  as  well,  but  is  obsured  by  the  dip  in  Et  due  to 
the  peak  in  ng. 

The  presence  of  a  significant  electron  pressure  gradient  effect  makes  it  impossible  to 
infer  density  with  the  present  probe  data  alone.  Such  an  inference  would  be  possible  under 
the  assumption  of  an  isothermal  plasma,  in  which  case  Equation  (4)  woidd  become,  for 
each  spatial  point,  a  differential  equation  in  ne(z).  The  actual  solution  of  the  whole  set 
of  these  would  be  a  cumbersome  iterative  process,  probably  not  justified  in  view  of  the 
uncertainty  of  the  isothermal  assumption. 

It  has  been  possible,  however,  to  obtain  an  informative  and  significant  distribution  of 
ne{z,i)  for  early  times,  specifically  during  the  breakdown  phfise.  This  is  shown  in  Figure 
17.  The  most  important  result  here  is  that  the  ionization  proceeds  unusually  slowly.  Full 
plasma  density  is  not  attained  for  2.5  w-seconds  and  during  this  time  the  discharge 
current  has  reached  half  of  its  peak  value.  Plasma  resistivity  2md  the  associated  diffusion 
of  the  field  are  high,  and  efficiency  is  accordingly  degraded. 

9  Simulation  of  Breakdown  Phase 

In  order  to  better  understand  the  physics  of  the  breakdown  of  the  neutral  gas  propellant  on 
the  face  of  the  PIT  coil,  a  computational  simulation  of  the  problem  has  been  undertaken. 
This  study  has  been  limited  in  scope  to  the  formation  and  heating  of  the  plasma,  and  the 
diffusion  of  magnetic  flux  through  it  during  ionization  and  current  sheet  formation. 

Any  plasma  simulation  that  includes  in  the  maximum  detail  all  the  physics  that  can 
plausibly  be  significant  in  a  problem  such  as  this  will  be  massive  in  terms  of  the  size  of 
the  code,  computational  facilities  required,  and  programmatic  effort.  A  few  such  codes 
exist.  However,  in  instances  where  the  area  of  interest  is  somewhat  restricted,  and  where 
physical  conditions  allow  reasonable  approximations  to  some  peu-euneters  to  be  made,  or 
others  to  be  neglected,  a  relatively  simple  program  cim  supply  significant  emswers. 

In  the  present  instance,  magnetohydrodynamics  has  been  neglected  altogether,  under 
the  assumption  that  plasma  motion  during  the  breakdown  phase  will  not  affect  the  process 
materially.  It  is  assumed  that  the  plasma  is  transparent  to  its  own  radiation,  so  that  radi¬ 
ation  transport  is  not  an  issue.  Close  to  the  coil,  over  most  of  its  area,  spatial  dependences 
of  variables  are  mainly  axial  (along  z),  so  a  one-dimensional  calctdation  is  adequate.  Other 
simplifying  approximations  zue  discussed  in  the  separate  subsections  below. 

0.1  Computational  Grid  and  Circuit 

The  one-dimensional  grid  consists  of  50  cells,  each  4  mm  long,  along  the  z  axis,  the  direction 
of  acceleration.  The  magnetic  field  B  is  defined  as  directed  along  y,  and  the  current  density 
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and  applied  electric  field  along  x.  This  configuration  may  be  visualized  as  a  parallel-plate 
cluumel  of  unit  plate  separation  (x),  infinite  width(y),  and  20  cm  length(z). 

Parasitic  circuit  inductance  is  provided  by  extending  the  acceleration  channel  toward 
negative  z  (where  plasma  is  excluded),  and  applying  a  fixed  driving  voltage  at  its  end. 


9.2  Ionization 


Ionization  is  assumed  to  be  due  to  direct  electron  impact  from  the  groxmd  state,  and  also 
from  the  first  excited  state.  For  the  latter  process,  it  can  be  assumed  that  excitation  from 
the  ground  state  is  the  rate-limiting  step,  and  that  ionization  occurs  promptly  thereafter. 
The  cross  sections  for  excitation  and  direct  ionization  are  taken  from  Druyvesteyn  and 
Penning  [4].  Each  rises  linearly  near  its  threshold  (10.8  ev  and  15.7  ev,  respectively),  and 
since  the  plasma  temperature  is  expected  to  be  no  more  than  about  3  or  4  ev,  only  the 
linear  portion  of  each  cross  section  near  its  threshold  will  contribute  significantly  to  the 
ionization  rate. 

The  total  ionization  rate  due  to  either  process  is 

Tl  =  ReHo  <  av  > 


where  n*  is  total  electron  density,  Uo  is  density  of  neutrals,  and  <  <tv  >  is  the  mean  of 
the  product  of  electron  velocity  and  cross  section,  teJcen  over  the  spectrum  of  electron 
velocities.  In  terms  of  the  original  neutral  density  n^,  =  rio  -  Hj. 

<  (TV  >  may  be  evaluated  by  starting  with  the  differential  ionization  rate: 


dn  =  no  (t{E)  n{E)  ^  dE 


where  <r  and  n  are  now  explicit  functions  of  electron  energy  E.  Electrons  are  Maxwell- 
distributed: 

n{E)  = 


so  that 


V^27rm7(lfcT)i 

where  k  is  the  Boltzmann  constant  and  T  is  the  electron  temperature  in  ’’K. 
The  linearly  rising  cross  section  near  a  threshold  Eo  is 


<t{E)  =  AiE  -  Eo) 


where  A  is  obtained  fix>m  previously  cited  experimental  data  [4].  This  finally  gives,  as  the 
total  ionization  rate  for  a  process. 


n  = 


r  ElB  -  Bo)  0-^1" dB 
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which  integrates  to 


Eo 


,-EolkT 


When  direct  and  excitation-ionization  rates  are  combined,  using  Eo=15.7  ev  and  Ao  = 
1.9  X  10“**  m*/ev  (direct  ionization)  and  .Ei=10.8  ev  and  Aj  —  0.3i4o  (excitation),  as  well 
as  the  approximation  £o.i  >>  the  result  is 


=  1.0  X  10“*^n*no  ^/T  (l  -I-  0.20  m-^s"*  (6) 

where  T  is  expressed  in  ev. 

Possible  radiative  ionization  has  been  neglected  in  this  initial  formulation  of  the  prob¬ 
lem.  It  would  be  a  considerable  complication  over  what  has  been  done  here,  but  could  be 
included  in  further  refinements  of  the  code  if  results  indicate  it  to  be  a  significant  omission. 


9.3  Heating  and  Losses 

Plasma  heating  is  due  only  to  resistivity;  both  electron-ion  and  electron-neutral  resistivity 
terms  are  included.  Losses  are  due  to  radiation  and  ionization  cooling. 

9.3.1  Resistivity 

Plasma  (electron-ion)  resistivity  has  the  Spitzer  dependence  on  temperature,  with  the 
coulomb  factor  In  A  equal  to  approximately  5: 

r?p  =  10-®r“^  n  -m  (7) 

Electron-neutral  resistivity  is  calculated  using  cross-sections  from  Rosa  [5]: 

Tj„  =  1.7  X  10~^\/T  Tiofrie  fl-m  (8) 


9.3.2  Radiation 

Radiation  losses  sue  calculated  using  an  approximation  due  to  Griem  [6],  which  is  that 
for  atoms  or  ions  whose  first  excited  states  are  not  unusually  low  in  energy,  the  principal 
radiation  loss  is  through  the  resonance  lines;  the  total  loss  is  faurly  well  approximated  by 
doubling  the  resonance  line  loss.  The  resonance  lines  of  both  AI  and  All  have  nearly  the 
same  energy,  which  allows  the  simple  approximation  that  all  of  the  original  neutral  density 
is  radiating  at  the  same  rate.  With  constants  evaluated  and  T  in  ev,  this  power  is 

P,  =  2.0  X  10“®‘n,n,  (9) 

Here,  the  energy  of  the  resonance  line,  Pj  =  11.8  ev,  appears  in  the  exponent,  ria  appears  in 
this  expression  to  provide  the  approximation  that  both  the  neutrals  and  ions  are  radiating 
the  same  li.  ?  energy. 
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9.3.3  Ionization  cooling 

It  is  assumed  that  each  ionization  removes  &om  the  plasma  an  energy  equal  to  the  ioniza¬ 
tion  potential  plus  the  temperature  of  the  newly  equilibrated  electron. 


9.4  Field  Diffusion,  Current,  and 

As  noted  above,  the  source  voltage  is  assumed  constant  during  the  ionization;  for  the 
unit  width  channel,  this  is  represented  by  a  constant  applied  electric  field  E,.  Parasitic 
inductance  between  this  source  and  the  plasma  is  represented  by  an  additional  length  Ip 
of  the  channel.  These  may  be  converted  to  actual  bank  voltage  and  parasitic  inductance 
by  multiplication  by  2irr  and  2irftorlp,  respectively,  where  r  is  the  radial  position  of  the 
plasma  on  the  coil  face. 

Through  Faraday’s  law,  E,  may  be  equated  to  the  time  derivative  of  the  flux  in  the 
channel: 

where  Bi  =  =  0)  is  uniform  over  Ip,  the  parasitic  segment  of  the  channel. 

For  the  discrete  computational  grid,  the  flux  integral  is  replaced  by  a  sxunmation  over 
the  plasma  cells.  The  equation  is  then  integrated,  and  a  fiurther  rearrangement  gives  am 
expression  for  updating  Bt : 


1  r  50  • 

5,  =  -  E,t-AzJ^B„ 

‘P  n=l 


The  distribution  of  B  in  the  plaisma  is  governed  by  the  diffusion  equation 


dB  _  Ti  d'^B 
dt  fio  dz"^ 


(11) 


which,  when  cast  into  finite  difference  form,  gives  the  increment  of  B  in  the  t  th  cell,  per 
timestep: 

ABi  =  -2- 

/^o 

where  tj  =  rjp  +  rjn,  the  sum  of  the  electron-ion  and  electron-neutral  resistivities. 

The  current  density  in  each  cell  is 


Bj+i  -H  Bj-i  —  2Bi 

(Azy 


r.  —  -L 

Ho  I  2Az 


(13) 
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9.5  Computation  Algorithm 

At  ezich  timestep: 

•  The  Bi  in  the  plasma  channel  are  summed  to  get  a  new  flux  integral. 

•  This  integral  is  differenced  with  that  from  the  previous  timestep  smd  divided  by  At 
to  obtain  Ex,  the  electric  field  at  the  coil  insulator. 

•  Bi  is  evaluated  with  Equation  (10). 

Then,  at  each  cell,  and  within  the  particular  timestep: 

•  »;(no,ne,r)  is  evaluated  through  Equations  (7)  and  (8). 

•  A  new  array  Bli  =  Bi  +  ABi  is  evaluated  through  Equation  (12)  (diffusion  equation) 
and  the  new  rj. 

•  n*  is  updated  through  Equation  (6). 

•  Te  is  updated  through  the  energy  balance  equation 


^  ~  _  ^15.7 +  r 


(14) 


•  The  B  array  is  replaced  by  the  BI  array,  time  is  ad%'anced,  zind  the  cycle  repeated. 


9.6  Stability  of  Computation 

The  diffusion  equation  (Equation  (12))  is  made  nonlinear  by  the  dependence  of  rj  on  field 
quantities,  so  it  might  be  anticipated  that  the  ordinary  stability  criteria  for  the  first  order 
explicit  integration  scheme  employed  here  might  not  be  adequate.  This  has  proven  to  be 
the  case.  For  a  constant  diffusivity  (rj/fto)  the  timestep  requirement  for  stability  is  that 

At  <  0.5^4^)^ 

For  these  calculations,  a  timestep  limit  for  stability  has  been  found,  but  it  is  approxi¬ 
mately  ten  times  smaller  than  that  given  above.  For  the  parameters  of  the  computations 
done  here,  a  timestep  of  1.0  x  10~’  sec  is  adequate. 
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9.7  Program  Inputs:  Initial  Conditions 

The  propellant  is  assumed  to  be  initially  deposited  against  the  coil  with  a  distribution 
decreasing  exponentially  in  2,  with  a  scale  thickness  of  6  cm.  Its  peak  density  is  a  variable 
input  parameter,  and  it  is  1%  ionized  at  t=0.  Its  initial  temperature  is  set  at  0.1  ev.  For 
the  Mkl  thruster,  the  peak  cold  propellant  pressure  is  lOOp,  and  for  MkIV,  600^*. 

JSt  and  Ip  are  also  assigned  as  program  inputs.  Since  the  experimental  data  against 
which  these  computations  ,are  compared  is  taken  at  a  radius  of  27  cm  on  MkIV  and  40 
cm  on  Mkl,  Vq  =  1.7£«  on  MkIV  and  2.5£,  on  Mkl,  while  the  actual  parasitic  inductance 
is  Lp  =  2.1  X  10“*/p  on  MkIV,  and  3.1  x  10"*/p  on  Mkl.  Lp  may  be  determined  for  each 
thruster  from  its  initial  dB/dt  or  alternatively,  from  the  time  of  its  first  current  peak,  when 
the  principal  inductance  in  the  circuit  is  still  the  parasitic  component.  For  Mkl,  this  gives 
Lp  =  130nff  {lp  =  .04m),  and  for  MkIV,Lp  =  Z20nH{lp  =  .15m). 

For  both  thruster  data  runs,  the  bsuik  voltages  corresponded  to  electric  source  fields 
Eg  of  approximately  10000  V/m. 

9.8  Results 

Since  plasma  djTiamics  are  not  included  in  these  computations,  the  comp^isons  made  here 
are  restricted  to  early  times.  All  computations  were  terminated  at  a  time  of  2.0  ps. 

For  the  rates  of  ionization  and  radiation  given  in  Equations  (6)  and  (9),  Table  2  com¬ 
pares  Mkl  and  MkIV  data  with  computations.  The  parameter  fion  is  fractional  ionization. 


Mkl  MkIV 


calc 

exper 

calc 

exper 

Br{2ns) 

.35 

.35 

.11 

.13 

EoiO) 

2.5  X  10’ 

2.5  ±  1  X  lO’ 

1.5  X  10’ 

1.3  ±  .3  X  10’ 

fion 

.80 

— 

.02 

.1  ±  .05 

Te 

3.9  ev 

— 

2.1  ev 

— 

Table  2:  Calculated  and  measured  parameters  for  Mkl  and  MkIV  thrusters.  Ionization 
rate  constant  =  1.0  x  10”*^. 

Agreement  of  B,  values  at  2ps  and  Ef  at  f  =  0  is  seen  to  be  excellent;  fractional 
ionization,  for  which  experimental  values  were  only  obtained  for  MkIV,  is  somewhat  low^. 

’This  pressure,  while  higher  than  what  would  ideally  produce  I,p  in  the  1500-2000  second  range,  was  that 
employed  for  the  complete  E,  B,  and  J  mapping  of  Run  25. 

’Since  the  computation  assumes  a  constant  source  voltage,  the  experimental  values  given  for  Br  are  those 
projected  to  2/is  from  the  earlier  slopes 
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In  order  to  compensate  for  the  neglect  of  radiative  ionization  and  ionization  from  excited 
states  other  than  the  first,  the  ionization  rate  constant  in  Equation  (9)  was  ajusted  upward 
from  its  original  value  of  1.0  x  Table  3  presents  values  for  an  increase  of  the  rate 

constant  by  a  factor  of  20. 


Mkl  MkIV 


calc 

exper 

calc 

exper 

Br(2fi3) 

.33 

.35 

.11 

.13 

Bg(0) 

2.5  X  10^ 

2.5±1  X  10^ 

1.6  X  10® 

1.3±.3x  103 

fion 

1.00 

— 

.07 

.1  ±  .05 

Te 

3.9  ev 

— 

1.5  ev 

— 

Table  3:  Calculated  and  measured  parameters  for  Mkl  and  MkIV  thrusters.  Ionization 
rate  constant  =  2.0  x  10“'^. 

Here,  there  is  fairly  good  agreement  between  calculated  *ind  measured  ionization  frac¬ 
tions.  The  Mkl  ionization  for  this  case  is  100%,  which  is  in  accord  with  the  obserration  in 
that  experiment  that  optical  spectra  from  neutrals  were  entirely  absent. 

The  different  electron  temperatures  and  conductivities  displayed  by  Mkl  and  MkIV 
result  in  differing  initial  current  sheet  thicknesses.  Br{z)  in  any  case  tends  toward  an 
exponential  ditribution;  the  clearest  display  of  diffences  in  diffusion  depth  is  then  a  semi¬ 
log  plot  for  which  the  slope  of  the  distribution  curve  is  proportional  to  the  current  layer 
thickness. 

Figure  18  shows  plots  of  log(Br)  vs.  z  for  the  two  thrusters.  The  current  sheet  thick¬ 
nesses  are  5.1  cm  for  Mkl,  and  9.0  cm  for  MkIV. 

10  Conclusions 

•  Prior  work  with  the  Mkl  PIT  in  the  ringing  discheurge  mode  has  demonstrated  the 

feasibility  of  the  concept,  with  fairly  good  efficiency  in  the  range  of  interest  for 
E^arth  orbit  missions.  . 

•  While  clamped  dischsirge  operation  is  expected  to  improve  thruster  efficiency,  this 
has  not  yet  been  demonstrated,  either  with  the  Mklll  or  MkIV  thrusters. 

•  The  present  analysis  of  Mkl  and  MkIV  first  quarter-cycle  data,  together  with  a  nu¬ 
merical  model  developed  to  describe  initial  plasma  breakdown  (ionization  and  current 
sheet  formation)  has  shown  the  crucial  importance  of  circuit  parasitic  inductance. 
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•  The  large  indicated  parasitic  inductance  of  MklV,  probably  due  mostly  to  internal 
capacitor  inductance  associated  with  diode  modifications,  is  not  an  essential  limita¬ 
tion  of  the  PIT  concept,  and  could  readily  be  brought  to  a  V2due  lower  than  that  of 
Mkl  in  a  new  design. 

•  It  has  been  understood  throughout  the  history  of  pulsed  thruster  research  that  par¬ 
asitic  inductance  is  detrimental  to  efiBciency,  through  the  theorem  that  eflBciency 
can  never  be  greater  than  the  fractional  increase  in  inductance  during  the  acceler¬ 
ation  stroke.  It  has  been  discovered  in  this  work  that  a  more  serious  consequence 
of  parasitic  inductance  is  its  effect  on  initial  plasma  breakdown  and  current  sheet 
formation. 

In  fact,  the  disastrous  loss  of  ionization  rate  incurred  through  high  parsisitic  induc¬ 
tance  is  sufficient  to  explzdn  the  poor  current  sheet  formation  and  low  thrust  efficiency 
of  all  thrusters  since  Mkl. 
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THRUSTER  EFFICIENCY  vs.  Up  FOR  DIFFERENT  THRUSTERS  (ARGON) 
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FIGURE  9.  MK  IV,  AZIMUTHAL  ELECTRIC  FIELD  vs.  TIME 
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FIGURE  11.  MK  IV, 
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FIGURE  16.  MK  IV, 
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